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Abstract 

In this paper, we present a new approach to the Kratky-Porod Model (KP) of semiflexible poly- 
mers. Our solution to the model is based on the definition of a generating function which we use to 
study the statistical mechanics of semiflexible polymer chains. Specifically, we derive mathematical 
expressions for the characteristic function, the polymer propagator and the mean square end-to-end 
distance from the generating function. These expressions are valid for polymers with any number 
of segments and degree of rigidity. Furthermore, they capture the limits of fully flexible and stiff 
polymers exactly. In between, a smooth and approximate crossover behavior is predicted. The 
most important contribution of this paper is the expression of the polymer propagator which is 
written in a very simple and insightful way. It is given in terms of the exact polymer propagator of 
the Random Flight Model multiplied by an exponential that takes into account the stiffness of the 
polymer backbone. 
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I. INTRODUCTION 



Since the ground-breaking ideas by Kratky and Porod in 1949, 1 theoretical studies of 
semiflexible polymers based on the KP model have been abundant. In this model, the poly- 
mer chain displays resistance to bending deformations. This resistance is modeled using a 
free energy that penalizes bending the polymer backbone. The free energy depends on pa- 
rameters (elastic constants) that are a consequence of many short-range monomer-monomer 
interactions. Explicitly, for the continuous version of the KP model 2 (called the Wormlike 
Chain Model, WCM) the free energy is 



where R (s) is the vectorial field that represents the polymer chain, s is the arc of length 
parameter, L is the contour length of the polymer and e is the bending modulus. In addition, 
the local inextensibility constraint |dR(s) /ds| = 1 must be satisfied. 

As a consequence of the bending rigidity, a semiflexible chain is characterized by a persis- 
tence length (proportional to the bending modulus) such that, if the length scale is shorter 
than the persistence length, then the chain behaves like a rod while, if the length scale is 
larger than the persistence length, then the chain is governed by the configurational entropy 
that favors the random-walk conformations. 

The local inextensibility constraint has not allowed researchers to find an exact solution 
to the model. Indeed, the constraint |dR(s) /ds| = 1 is written using a Dirac delta distri- 
bution in infinite dimensions. Depending on how the constraint is written, |dR(s) /ds| = 1 
or (dR(s) /ds) 2 = 1, we get an Edwards Hamiltonian that is non-analytic or non-linear, 
respectively. Consequently, there is no exact solution of this model at present. However, a 
few properties like the first few moments of the distribution of the end-to-end distance 2 ' 3 are 
known exactly. 

Indeed, many researchers have addressed the Kratky-Porod and other models of semi- 
flexible polymers with the purpose of understanding the statistical behavior of this kind of 
polymers. Among the many theoretical treatments of semiflexible polymers, let us start 
by mentioning two very important contributions done in the 1950s: the work by Daniels 4 
who developed expansions of the polymer propagator for a semiflexible polymer in inverse 
powers of the number of segments and the classic paper by Benoit and Doty 5 who obtained 
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the exact expression for the average end-to-end distance squared and radius of gyration. 
During the following two decades, the field of statistical mechanics of semiflexible polymers 
saw a substantial growth thanks to the seminal contributions of many researchers. For ex- 
ample, Fixman and Kovacs developed a modified Gaussian model for stiff polymer chains 
under an external field (external force). 6 In this approach, they computed an approximate 
distribution for the bond vectors from which they were able to compute the partition func- 
tion and average end-to-end vector. An alternative approach was proposed by Harris and 
Hearst 7 who developed a distribution for the continuous model from which they were able 
to compute the two-point correlation function and, consequently, the mean-square end-to- 
end distance and radius of gyration. A reformulation of the KP model using field-theoretic 
methods was developed by Saito and coworkers 2 who computed exactly different averages 
of the end-to-end distance and tangent-tangent correlation function. In addition, Gobush 
and coworkers 8 developed an asymptotic expression for the polymer propagator in inverse 
powers of the number of segments and, Yamakawa and Stockmayer 9 addressed the first order 
corrections to the end-to-end distance squared and second virial coefficient due to excluded 
volume interactions. Semiflexible polymers were also extensively studied by Freed 10 who 
employed field-theoretic methods to study their statistical behavior. A few years later, he 
developed a modified Gaussian distribution approximation to the continuous version of the 
KP model 11 which has been re-derived using different mathematical methods by Lagowski 
and coworkers, 12 and Winkler and coworkers. 13 Similar results were obtained by Zhao and 
coworkers 14 who also studied the effect of an external field. Field-theoretic methods have 
also been used by Bhattacharjee and Muthukumar 15 who employed the Edwards-Singh self- 
consistent approach to obtain an effective Gaussian Hamiltonian from which they computed 
the average end-to-end distance squared. 

In recent years, the advent of new experimental methods that have allowed researchers to 
manipulate single molecules has generated new momentum in the area of statistical mechan- 
ics of semiflexible polymers. 16 Indeed, new analytical approaches to semiflexible polymers 
have been developed by Marko and Siggia, and Kroy and Frey 17 who derived the force ver- 
sus elongation behavior predicted by the WCM, Hansen and Podgornik 18 who developed a 
non-perturbative 1/d expansion (d being the dimension of the embedding space), Whilhem 
and Frey 19 who computed the polymer propagator for polymers with large bending rigidities 
and Winkler 20 who computed the same quantity for any value of the stiffness of the poly- 
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mer backbone using the Maximum Entropy Principle. Finally, it is worth mentioning one 
alternative approach to semiflexible polymers (Dirac chains) developed by Kholodenko. 21 

The aforementioned extensive list of approximate treatments of the KP model, though 
not exhaustive, does show two major trends. On the one hand, many studies computed 
different statistical properties as perturbative corrections to the rod-like and flexible limits. 
One example of this statement is the seminal paper by Gobush etal. where the polymer 
propagator was computed as a perturbation expansion with respect to the Gaussian chain 
model (see Eqs. (22) and (23) in Ref. 8). On the other hand, many researchers have targeted 
a few physical properties of the KP model and optimized the approximations to capture these 
properties exactly or very accurately. Two examples of this statement are: the paper by 
Bawendi and Freed, 11 who optimized their treatment of the KP model (continuous version) 
to reproduce the bond-bond projection expectation value exactly and obtained a modified 
Gaussian distribution for the end-to-end distance (see Eq. (13) in Ref. 11), and the paper 
by Winkler etal. 13 who used a very fundamental concept, the Maximum Entropy Principle, 
to compute different ensemble averages and correlation functions (the constraints used to 
optimize the use of the Maximum Entropy Principle are given by Eqs. (3.25) to (3.28) in Ref. 
13). These observations lead naturally to the following two questions . Firstly, we should 
ask "why do we need so many different approximations to solve the same model?" and, 
secondly, "can we find one approximation capable of predicting all the statistical properties 
of the model in such a way that corrections to these results can be computed in a systematic 
and perturbative fashion ?". Observe that the answer to the latter is affirmative for the 
problem of a single flexible polymer with excluded volume interaction. 22 

In this paper we address the second question. In other words, we address the KP model 
with the purpose of finding one "ground" state capable of reproducing all the statistical 
properties of the model approximately. Furthermore, we require that our "ground" state 
captures all the statistical properties of the rod-like and flexible limits exactly and respects 
the local, not global, inextensibility constraint. In the crossover region, the solution is not 
exact but, we require that it displays the correct physical features as described by other 
treatments of the model and, moreover, it should be amenable to systematic and controlled 
corrections calculated in a perturbative manner around the "ground" state. In this first 
paper, we set the Fixman parameter to zero and leave the excluded volume problem for a 
future publication. 
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The most important prediction of this calculation is the polymer propagator which turns 
out to be proportional to the polymer propagator of the Random Flight Model 23 multi- 
plied by an exponential arising from the stiffness of the polymer backbone described by the 
parameter a (see text for details), 



where n is the number of segments and R is the end-to-end distance. The mathemati- 
cal simplicity of this expression makes it a very good candidate for the development of a 
perturbative treatment of the KP model. 

This paper is organized as follows. In section II, we propose a generating function and 
obtain an approximate analytical expression for it. This section also contains the only 
approximation of our calculation. Afterward, we use the generating function to calculate 
the characteristic function, mean squared end-to-end distance and polymer propagator of 
the model. In Section III we discuss the results of our calculations which are valid for any 
value of the stiffness of the polymer backbone and length of the polymer chain. For the 
purpose of making our presentation more balanced and objective, we compare our results 
with those obtained by other researchers. Section IV contains the conclusions of our work. 
The details of some mathematical calculations are presented in the appendix. 

II. THEORY 

A. Review of the Kratky-Porod Model and Evaluation of the new Generating 
Function 

Let us consider a polymer chain as a set of n bond vectors (u , u 1; u n _i) connected in 
a sequential manner. Furthermore, let us assume that the length of each bond vector is 4 
(=Kuhn length) and that pairs of consecutive bond vectors try to be parallel to each other. 
This orientational interaction is modeled with a Boltzmann weight given by the following 
expression 6 ' 15 



P (R) ~ P RFM (R) exp [R 2 /2an 2 ] , 
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where e and fcgT are the bending modulus and thermal energy, respectively. In addition, 
we take into account the local inextensibility constraint with the following term 2 



Equations Q and Q define the KP model completely and all the statistical properties 
of the model such as the characteristic function, single chain structure factor and other cor- 
relation functions, probability distributions like the polymer propagator and their moments 
can be calculated. The evaluation of these statistical properties can be easily done using the 
following generating function 



This definition shows that ^kj P is an auxiliary tensor. It represents a standard external 
source consisting of two terms: the first one, (Jfc), represents a dipolar coupling with the 
k-th bond vector and, the second one, (q), is also a dipolar coupling that is non-zero only 
when the k-th bond vector is between the j-th and p-th ones. The 0(z) function is the 
Heaviside step function. 24 

If the generating function given by Eq.(PJ) is known, then all the aforementioned 
statistical properties can be obtained from it by differentiation with respect to Jfc 
and/or by assigning specific values to Jfc . For example, the characteristic function is 
C ({iq ■ [Q(k — j) — Q(k — p)]}) from which the polymer propagator is computed as the 
inverse Fourier transform of C({zq}), the tangent-tangent correlation function, (uj • u&), 
is computed as the second derivative with respect to Jj and Jfc which has to be evalu- 
ated at J = Ji = ....J n _i = 0. Other statistical quantities can also be computed easily. 




(3) 





where ^kj P is 



y kjp = j fc + zq • [e(k - ]) - e(k - P )} . 



(5) 
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Consequently, our first step is to evaluate C ({^kj P }) whose explicit form is 



{n-2 n-1 ^ 

k B k=0 k=0 ) 

(6) 

where z is defined such that C ({0 k j P }) = 1- After writing the delta distributions in Eq.® in 
terms of their Fourier representations, 25 a Hubbard-Stratanovich transformation 26 is needed 
in order to avoid dealing with the quadratic term in arising from the exponential repre- 
sentation of the delta distributions. As a result, we get 



C ({<W) = -J 



n-l 

n 

i=0 



duid 




C ({*«.}) 



n— 1 



| J duidXid^i 



,i=0 



cxp 



n— 1 



n—2 g n—1 



k=0 



lnX k + i-X k 



E- - E^ 

fc=0 v fc 7 fc=o ft . 



(7) 



where and are dummy variables arising from the Fourier representation of the delta 
and the Hubbard-Stratanovich transformation, respectively. Observe that after redefining 
^kjph — > ^kj P and u k /l k —* u k , the dependence of Eq.JZJ) on the Kuhn length goes away. 

We now proceed to compute the integrals over X k and u^. The A^-integrals are straight- 
forward, the result is 



n-1 



i=0 



exp 



" -1 *l 3 



n-1 



Ell _ _ V 
A- o. 2^ 



k=0 



k=0 



In X k + i-X k 



n-1 



where a constant prefactor independent of $j has been neglected because of the definition 
of z in Eq.((ni). 

The u fc -integrals are exactly doable using saddle point which is given by the following 
expression 



k B T 



U k = Uq - 



fc-1 



mjp 



i$ m ) (k - m) 



(9) 
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where the following constrain must be fulfilled 



n-l 



(10) 



m=0 



Replacing Eqs.ljHJ) and (jHJ into Eq.(j2j), and expressing the constraint as a delta distribution 
in the exponential representation, the generating function takes the form 



C ({<W) 



nn 
1> 



kjp) 
n—1 n—1 



du 



n—1 



i=o 



cxp 



n—1 



i=o 

n-l 



sin(|$ 3 -|) 
(l*il) 



fe=0 s=0 s=0 

where the kernel K^ s and the parameter k are defined as 



(11) 



1 _k_e^ {s _ k) 

n n 



and k 



k B T 



(12) 



We refer the reader to Appendix A for some mathematical details of this part of the calcu- 
lation. 

The last step involves the evaluation of the ^-integrals. Observe that this part of the 
calculation is not exactly doable thus, we approximate it as follows. First, we assume that, 
for a fixed value of s or k, the element of Kk, s that contributes the most to the integral is 
the diagonal one. Thus, we approximate Kk, s by its diagonal elements. The contributions 
arising from non-diagonal elements will be treated in a future publication. We also perform 
a pre-averaging approximation. Consequently, the kernel s becomes 



K, 



k, s 



k . . 

1 /hs 

n , 



(13) 



This is the only approximation of our work. In the next section, we will show that this 
approximation captures the flexible (Random Flight Model) and rod-like limits exactly when 
k — > oo and k — > 0, respectively. This approximation also provides a very good crossover 
behavior for intermediate values of the parameter k. 

Equation ([TTfy transforms the integrand of the u-integral into 



n-l 

n 

8=1 



sin(|$J) 
d$ s , , ; exp 



-iu • ($ s + i^ sjp ) - ^($ s + iV sjP Y 



(14) 
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where, for the purpose of simplicity, we have introduced a new parameter a = | (l + i) k. 

Finally, replacing Eq. lfT^) in Eq. lfTTJ) we obtain the most general expression of the gener- 
ating function C ({^kj P })- In the next section we evaluate some physical properties of the 
KP model using this generating function. 



B. Characteristic Function and the Mean Square End-to-End Distance 

Let us start by computing the characteristic function. For this purpose, we make \1/ = zq 
in Eq.tJTTJ where q is the scattering wave vector. Then, the characteristic function and the 
mean square end-to-end distance can be calculated. The ^-integration takes the form 



] [ / d$ / d(cos{9)) / ci^sm($)$exp - — ($ 2 + g 2 ) + 
_ n Jo J -i Jo L 2 



3=0 



(15) 



+?q ■ u + anq • $ — iu ■ $] , 

whose calculation is straightforward. Choosing 9 to be the angle between the vectors 
u and $, and 9 q the angle between q and <3>, then the solution to the (^-integral is 
2ttIq (an$g |sin (9) sin (9 q )\) where Iq (x) is the Bessel function of second class. 27 Hence, the 
angular integral, Iq, reads 



I n = 2n I d (cos (6>)) exp [(— m$ + an§q cos (9 q )) cos (9)] I (an§q |sin (9) sin (9 q )\) 



sin ( $y (u + ianq) 2 ) (16) 

= 47T . 

<&y (u + ianq) 2 

When we replace this expression into Eq. (fT5^l we obtain the following form for 7$ 



sin I $\/(u + ianq^ 2 



Z" 00 V v / / an \ 

Iq,=4ir c/$sin($) v exp — — ($ 2 + q 2 ) + iq • u . (17) 

• ' (u + «anq) 
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The solution of this integral is 



- sinh \Jk\/ (u + icmq) 2 ^ 
(u + zcmq) z 



\2 



Air J v ; = ^exp ^ -^+zqu q 2 . 18 

A/ 2cm H ^ M 2an 2cm H 2 H ' v J 



Consequently, the characteristic function can be written as follows 



47T 



7T 



2 cm 



OO />1 

diz / d(cos(0)) 
u 7-i 



sinh (u + ianq) 



u + icmq) z 



(19) 



where A/" is the norm which is defined in such a way that C (0) = 1. 

In the limit of a — > oo, Eq. ([TT?l) approaches the asymptotic limit given by the formula 



sin (q) 



q 



(20) 



which is the characteristic function of the Random Flight Model. 28 
On the other hand, if a — ► 0, 



Cq^o(q') ~ / du I dVtu 2 n exp( — — — h iqncos (6) ) sinh n (u/an) 
Iq J V 2a 

sin (gn) 



(21) 



and the Rod-like limit is recovered. 29 

An expression for the entire range of values of a is obtained from Eq. (JTHj)- The first step 
is to compute the integral over the angles, that is 



C Q (q) = 2tt / d(cos(0)) 



sinh ( — 



u + zcmq) 



(u + ianqY 



(22) 



A variable transformation v 2 = (u 2 — (anq) 2 + 2iuanq cos (9)) allows us to rewrite Cn(q) 
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as follows 



(u-\-ianq) 

Cn(? ) = ^( aJi )-/^- d „E^M. (23 ) 
iqu J v n 1 

an 

Replacing the hyperbolic sine by exponentials and, after expanding the resulting binomial, 
we can rewrite Ca(q) as follows 

which is easily calculated in terms of the Incomplete Gamma function. 27 The result is 



Cn(q) 



2ttz 



(2an) 1 - n J2(- 1 T\ " \ (2k-n) n ' z Im^T [2- ni (2k-n) 

1 U l—n \ k 



k=0 



(u + ianq)\ 



an 



(25) 

Going back to C(q), we note that the range of integration can be expanded to all real 
values of u by adding C(q) to its complex conjugate to get an integrand invariant under the 
transformation u — > —u. Then, we integrate it by parts to get 



(3n/2+l) / 2 \ (™/ 2+1 ) " / n 



Im 



du [u + iqan) 1 n exp 



u 

~2a~ + 



n-2k 



an 



u + iq(n — 2k) 



(26) 



The calculation of the integral in Eq. ([2~K|) is straightforward. The result is 



C(q) 



3(n + l) 

47T 2 



Re ( exp 




iq (n — 2k) 



nxn 



(27) 

where U(a, b, c) is Kummer's function 27 which can be rewritten in terms of Parabolic Cylin- 
der, Whittaker or other functions. Its mathematical properties like recurrence relations and 
integral representations are very well known and tabulated. 27 But, although its mathemat- 
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ical properties are many, we were unable to compute the real part in Eq. lJzTJ) in a compact 
mathematical form. 

Kummer's function has a cut on the negative real axis. Consequently, the evaluation of 
the norm N and the mean square end-to-end distance (R 2 ) cannot be obtained from the 
Taylor expansion of Eq. f2*T|) in powers of the wave vector q. Taking into account that both 
real and imaginary parts of U(a,b, — \c\) have well-known expressions 30 we used Eq. (|26|) . 
Firstly, we interchanged the sum with the imaginary operator and carried out the sum 
which gives a hyperbolic sine to the nth power. Afterward, the integrand was expanded in 
powers of the wave vector to third order and the integrals evaluated. The results are 

(27i) 3n/2+1 

N = (28) 

(an) 

and 



(2ir) 3n/2+1 

;R 2 > = -77 \ ?M { {n 2 ~ n) I n n Z 2 2 - n 2 [2a (n - 1) - 1] J£_ a - (n - 1) (n - 2) } 



N (an) 3n/2 - 3 



where 



(29) 



cxp 



an 
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^JHn m a 



2P+m/2-l/2 



cos 



7rp\ / V 



im/2 



2 M p/2 / r (™/2 + 1/2) 



b/2] 

+2 exp 



fc=0 



(p - 2kf 
2an 2 



P 
k 



2'2 7 2an 2 



C. Evaluation of the Polymer Propagator 



We now proceed to evaluate the polymer propagator. For this purpose, the starting 
point is the Fourier transform of the characteristic function given by Eq. lf^ . Firstly, we 
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interchange the order of integration as follows 



/+oo 
duexp 
-oo 



u 2 n — 2k 
2a an 




P(R)=R- 1 N~ 1 -n^ 2 -V (—) 
K ' 2 \anj 

Im |y dq sin (qR) (u + iqan) 1 ^ 11 exp [iq (n — 2k)] | , 



(30) 

and evaluate the imaginary part of the g-integral assuming that n is even. After replacing 
sin (qR) by its definition in terms of exponentials, the propagator can be written in terms 
of the following integrals 

2(an) l^Jo (q — iujan) ) 

which can be solved if we consider q to be a complex variable and apply Cauchy's theorem. 

The evaluation of the VF-integrals requires special attention to the location of the pole 
q = iu/an because u can be positive or negative. For reasons of notation, we define 
W*, W^, Wi, W2 such that the subindices 1 and 2 indicate a solution valid for u > 
and u < 0, respectively. 

We start by evaluating W + using Cauchy's theorem. If (n — 2k + R) > 0, we choose the 
contour of integration CI given in Fig. 1. The result for u > is 



Wt (R) = - , * 1( {(n - 2k + i?)- 2 exp 

2 (an) (n — 2)1 { 

On the other hand, if u < 0, then (R) = 0. 

The other possibility is (n — 2k + R) < 0. Then, the contour of integration in Fig. 1 
should be C2. This implies that (R) = and 



an 



(n-2k + R) 



(32) 



W+ (R) 



7Y 



2(an) n ' 1 (n-2)\ 



[n — 2k + R) n 2 exp 



an 



(n-2k + R) 



(33) 



The evaluation of W is done in a similar way. If (n — 2k — R) > 0, then the contour 
of integration should be CI and the results are (R) = and W{(R) = W^(—R). 
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Otherwise, if (n — 2k — R) < 0, then the contour is C2 and the results are W 1 (R) = and 
W^{R) = -W+(-R). 

Substituting these results into Eq. (|30j) and evaluating the following integrals 



and 



A k (R) = (n-2k + R) 



n-2 



dtexp 



in 



2k + R) n ~ 2 ^na/2 exp [R 2 /2an 2 ] 



Rt/aLn 

2a 1 

l-erf[ -i- 1 



B k (R) = (n-2k- R) n ~ 2 j dt exp 

J — oo 

= (n - 2k - R) n ~ 2 v / vra/2exp [R 2 /2an 2 ] 
then, the polymer propagator given by Eq. (p?U| becomes 



2an J 



h Rt/an 

2a 

1 - erf ( S-) 



2an J 



(34) 



(35) 



7r (3n/2+l/2) a 3(l-n)/2 2 n/2+l/2 

P {R) = nW(n-2)\N 6XP 



R 2 



2 -. [(n-fl)/2] 



2an 2 



E 

fc=0 



n 



(n - 2Jfe - i2) 



n-2 



(36) 



This expression, which was obtained assuming that n was even, is easily extended to any 
number of segments by analytic continuation. 

A comparison of this result with the one of the Random Flight Model 23 shows that we 
can write the final expression for the propagator in a very simple and insightful way 



P(R) 



a 3/2 (27r) (3n/2 + 3/2) 



N {anf n/2 

where the RMF stands for Random Flight Model 



Prfm (R) exp 



R 2 



2an 2 



(37) 



III. RESULTS AND DISCUSSION 

The polymer propagator described by Eq. ([37]) clearly shows that our treatment of the 
KP model captures the limit of fully flexible chains (a — > oo) exactly. Indeed, in this limit, 
Eq. (|37]) approaches the well-established expression of the Random Flight Model. Further- 
more, Eq. ([37]) clearly shows how the stiffness of the polymer backbone modifies the propa- 
gator of the RFM. In the other limit (a — > 0), the exponential grows. Consequently, those 
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configurations of the polymer chain with large end-to-end distance have higher probability 
of been realized than the configurations with small end-to-end distance. This is the correct 
physical behavior and is a consequence of the free energy penalty related to the formation 
of hairpins. But, although the exponential function grows in this limit, the propagator of 
the RFM puts an upper limit to the possible values of the end-to-end distance. This limit is 
the number of segments of the chain. This is clearly shown in the upper limit of the sum in 
Eq. (|36|) . Thus, the polymer chain satisfies the local inextensibility constraint. Another very 
important consequence of Eq. ()37|) is its mathematical simplicity which makes this expression 
of the polymer propagator a very good starting point for a perturbative treatment of the 
KP model as mentioned before. 

We plot the normalized radial distribution as a function of the end-to-end distance in Figs. 
2, 3, 4 and 5. The values of the parameters chosen are: n =6, 10, 20 and 30 Kuhn segments, 
and a =0.01, 0.02, 0.03, 0.05, 0.1 and 0.75. The figures clearly show that the location of the 
peak in the polymer propagator (multiplied by R 2 ) moves toward larger values of R when the 
stiffness of the polymer backbone increases. This behavior is in good qualitative agreement 
with previous results arising from computer simulation studies 19 and theoretical approaches 
based on the Maximum Entropy Principle. 20 This is the correct result because the stiffer the 
polymer backbone, the higher the energetic penalty to bend the chain. Consequently, those 
configurations of the macromolecule with small end-to-end distance will be more and more 
hindered as the stiffness increases, while those configurations with large end-to-end distance 
should be more and more favored. Therefore, the peak should shift toward larger values of 
R when the stiffness increases as shown by the figures. 

Figures 6 and 7 show the polymer propagator for polymer chains with 6, 10, 20 and 30 
Kuhn segments and a fixed value of the semiflexibility parameter. Theses figures show that 
the longer the polymer is, the more it behaves like a flexible chain since the location of peak 
(=end-to-end distance divided by the number of segments) moves toward smaller values. 
In other words, the longer the polymer is, the less relevant the stiffness of the backbone 
becomes. 

The two aforementioned effects can be easily rationalized in terms of Eq. (|37]) . Observe 
that the stiffness parameter a and the length of the chain n do not appear as independent 
variables but, as the product an 2 . Therefore, making the chain longer is equivalent to 
increasing the semiflexibility parameter thus making the polymer more flexible. 
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Another interesting result is the limit of very long chains. It is known that in this limit 
the polymer propagator of the RFM approaches the one of the Gaussian Chain Model. 28 
Thus, Eq. (jHZ| should approach the following expression 



P (R) ~ exp 



3R 2 R 2 1 
+ 



2n 2an 2 



exp 



3R 2 
' 2n 



3an 



exp 



3R 2 



2n(l 



K) 



(38) 



where we have defined an "effective" Kuhn length l e K (in units of the bare one, Ik) as follows 



i 



(39) 



( ^ 3an ) 

Observe that this renormalized Kuhn length approaches the bare one in the limits of very 
long , n — > oo, or very flexible chains, a — > oo. This is the expected result and can be 
rationalized using the arguments presented in the previous paragraph. In addition, note 
that as the chain becomes shorter or stiffer, Eq. (f3!?|) shows that the effective Kuhn length 
grows. Therefore, the effects of semiflexibility become more relevant. 

Figures 8 and 9 show the behavior of C (q), Eq. fjTTj) . as a function of the wave vector q 
for different values of a and two values of n (10 and 20). The figures clearly show that our 
treatment of the KP model predicts a smooth transition from the rod-like to the flexible 
limit. Moreover, our results capture the continuous change in the qualitative behavior of 
this function which changes from a monotonically decreasing function in the flexible limit to 
an oscillating function in the rod-like regime. In addition, our computations predict that the 
decrease of the characteristic function for small values of q should be faster for stiff polymers 
than for flexible ones. This is a consequence of the fact that rigid polymers have a larger 
mean squared end-to-end distance than flexible ones for a fixed chain length. 

The effect of chain length on the characteristic function is shown in Fig. 10 for the 
particular cases of chains with 6, 10, 20 and 30 Kuhn segments, and a fixed value of the 
semiflexibility parameter a(— 0.01). This figure shows that the longer the polymer is, the 
more it behaves like a flexible chain since the characteristic function approaches the one of 
a flexible polymer chain. 

Figures 11 and 12 show the mean squared end-to-end distance (R 2 ) as function of the 
semiflexibility parameter a for different chain length. Specifically, Fig. 11 compares our 
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approximate solution, Eq. (f29j) . for n = 8 with the exact one. 5 This figure shows that our 
approximate solution captures the limits of flexible and rigid polymers exactly and provides 
an very good behavior in the crossover. Figure 12 shows the behavior of (R 2 ) as a function 
of a for chains with 5, 8, and 10 Kuhn segments. 

Finally, we compute the exponent 2v in the relationship (R 2 ) = n 2u . Figure 13 shows 
the plot of /n{(R 2 )} as function of the number of segments, ln(n), for two values of the 
semiflexibility parameter, a = 0.001 (rigid) and a = 0.75 (flexible). The numerical values of 
the slopes are 1.95 and 0.98, respectively. These results are in excellent agreement with the 
expected values for the exponents of rodlike and fully flexible chains. Thus, our approach also 
captures the fractal dimensions of both limiting behaviors correctly. We have also studied 
intermediate values of the semiflexibility parameter and observed a smooth crossover from 
the rod-like to the flexible limit. 

IV. CONCLUSIONS 

In this paper we have proposed a new treatment of the KP model based on a generating 
function. The advantage of our approach is two fold. Firstly, the evaluation of the most 
relevant statistical properties of the model is straight forward since they can be obtained 
as derivatives or integrals of the generating function. This makes the evaluation of the 
generating function the crucial step for the solution of the KP model. In our treatment of 
this model, we were able to devise one approximation that was able to capture the flexible and 
rigid limits exactly and, moreover, was able to provide a smooth crossover behavior between 
the two aforementioned limiting regimes. This crossover has all the correct qualitative 
features as discussed in the previous section. Furthermore, this approximation respects the 
local, not global, inextensibility constraint. Secondly, our treatment of the KP model was 
able to provide a good, though not unique, "ground" state around which a perturbative 
treatment can be developed. We speculate that this perturbative calculation should be able 
to correct the consequences of our approximation in a systematic and controlled fashion. 
In other words, this perturbative analysis should account for the corrections arising from 
the neglected non-local terms in the kernel and deviations of the diagonal terms from the 
pre-averaging approximation used. 

The results presented in the paper are valid for polymers with any number of segments and 
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value of the stiffness of the polymer backbone. Specifically, we computed the characteristic 
function, the mean squared end-to-end distance and the polymer propagator of the KP 
model. Other averages like the radius of gyration and/or correlation functions like the 
single chain structure factor or tangent-tangent correlation functions can be computed in a 
straightforward manner, too. 

It is worth mentioning that the behaviors predicted for the characteristic function, poly- 
mer propagator and mean squared end-to-end distance are in very good qualitative and 
quantitative agreement with results arising from other studies of the KP and Worm-like 
Chain Model. This agreement makes us speculate that our future evaluation of the struc- 
ture factor, the radius of gyration and other correlation functions calculated within the 
present level of approximation should give good results at least at the qualitative level. 

Perhaps, the most important contribution of our work is the derived polymer propa- 
gator. The mathematical expression obtained for the propagator is very simple, compact 
and insightful. It becomes the exact expression of the Random Flight Model in the flex- 
ible limit. But, as the stiffness of the polymer backbone increases, the propagator of the 
RFM is modified by an exponential which depends on the semiflexibility parameter. As a 
consequence, the peak of the radial distribution function shifts toward larger values of the 
end-to-end distance when the polymer becomes stiffer. In the limit of infinitely stiff poly- 
mers, the propagator becomes a delta function centered at the total contour length of the 
chain. This is the correct result because the stiffer the polymer backbone, the higher the 
energetic penalty to bend the chain. This shift of the peak is accompanied by an increase 
in its height and a decrease in the width of the distribution. 

V. ACKNOWLEDGMENTS 

We acknowledge the National Science Foundation, Grant # CHE-0132278 (CAREER), 
the Ohio Board of Regents Action Fund, Proposal # R566 and The University of Akron for 
financial support. 



18 



Appendix A: SOME MATHEMATICAL ASPECTS OF THE EVALUATION OF 
THE GENERATING FUNCTION 



In this appendix we explain how the final expression of the generating function, Eq. (jTT|l . 
was obtained. For the purpose of clarity we define 

T k = ($ fc + M kjp ) , (Al) 

then, the w^-integrals in Eq.(j2J) can be evaluated exactly using saddle point which is the 
solution to the following equations 

Ul - u = -kT , (A2) 
u„-i - u n _ 2 = KT n _i, (A3) 

u fc+1 - 2u k + u fc _i = -nT k , fc = l..n-2, (A4) 

where k was defined in Eq.([T2}. The solutions to these equations are given by Eq.lJHl). More- 
over, the constraint for the sources, Eq. lfTUl) . is a consequence of these equations. Therefore, 
if this constraint is not satisfied, then Eq. © is not the saddle point solution. In addition, 
those terms appearing in Eq.Q can be rewritten as follows 

1 n-2 n-2 / k \ 2 



and also 



2k ^ k ^ ' 2 

fc=0 A;=0 \m=0 



n— 1 n—1 n—1 

K 



^u fe ■ T fc = -- ^ ^ T m |A; - m| T fc . (A6) 

fc=0 k=0 m=0 



We can simplify the sum of these two terms further if we note that 

n—1 n—1 n— 2 / k 

x m \k - m\ x fc + y^ r r 



fc=0 m=0 k=0 \m=0 
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n—1 n—1 n—2 fc n—2 n—1 

= ^^T m |A;-m|r fc + ^^r m r fe (72-l-£;)+^ ^ T m T fe (n-l-m) (A7) 

fc=0 m=0 fc=0 m=0 fc=0 m=k+l 

n—1 n—1 n—2 n—1 n—2 n—1 

= ^2^2r m \k-m\r k -^2 ^2 T m(m-k)T k + (n-l)^2^2T m T k 

fc=0 m=0 fc=0 m=fc+l fc=0 m=0 

where the last double sum vanishes because of the constraint, Eq. (fl"0|) . After changing the 
order of summation we arrive at 

n—1 n—1 n—1 k—1 

^ ^ Ym \k - m\ T fc - ^ T m (k - m) T fc . (A8) 

k=0 m=0 k=l m=0 

The first double sum, apart from the term k = 0, can be separated into two terms. One 
for m < k — 1 and another one for m > k such that the former cancels the second double 
sum in the previous expression. As a result we obtain 

n—1 n—1 n—1 

22 Y kT k + ^ ^ T m \k - m\ Tfc 

k=0 k=l m=k ( A q\ 

n-1 n-1 

= ^2 T ™( m ~~ k)e(m - k)T k . 

fc=0 m=0 

If we now use the constraint for the T fc 's, Eq. (fl"0|) . we can readily get the quadratic term 
in Eq.flTTl). 
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FIG. 2. Normalized polymer propagator 4nR 2 P (R) versus R/n for n — 6. Continuous line 
(a = 0.01) , dashed line (a = 0.02), dotted line (a = 0.03), dashed-dotted line (a = 0.05), 
dashed-dotted-dotted line (a = 0.1) and dashed-dashed-dotted line (a = 0.75). 



Comment: Figure 2, First Author: Marcelo Marucho, JCP 
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FIG. 3. Normalized polymer propagator 4nR 2 P (R) versus R/n for n = 10. Continuous 

line (a = 0.01) , dashed line (a = 0.02), dotted line (a = 0.03), dashed-dotted line 

(a = 0.05), dashed-dotted-dotted line (a = 0.1) and dashed-dashed-dotted line (a = 0.75). 
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FIG. 4. Normalized polymer propagator 4nR 2 P (R) versus R/n for n = 20. Continuous 

line (a = 0.01) , dashed line (a = 0.02), dotted line (a = 0.03), dashed-dotted line 

(a = 0.05), dashed-dotted-dotted line (a = 0.1) and dashed-dashed-dotted line (a = 0.75). 
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FIG. 5. Normalized polymer propagator 4nR 2 P (R) versus R/n for n = 30. Continuous 

line (a = 0.01) , dashed line (a = 0.02), dotted line (a = 0.03), dashed-dotted line 

(a = 0.05), dashed-dotted-dotted line (a = 0.1) and dashed-dashed-dotted line (a = 0.75). 



Comment: Figure 5, First Author: Marcelo Marucho, JCP 
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FIG. 6. Normalized polymer propagator 4nR 2 P (R) versus R/n for a = 0.01. Continuous 
line (n — 6), dotted line (n = 10), dashed line (n = 20) and dashed-dashed-dotted line 
(n = 30). 



Comment: Figure 6, First Author: Marcelo Marucho, JCP 
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FIG. 7. Normalized polymer propagator 4nR 2 P (R) versus R/n for a = 0.1. Continuous 
line (n — 6), dotted line (n = 10), dashed line (n = 20) and dashed-dashed-dotted line 
(n = 30). 
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FIG. 8. Characteristic function C (q) versus wave vector q for n = 10. 
Dashed-dotted-dotted line (a = 0) (the exact solution of rigid Model), 
dashed-dashed-dotted line (a = 0.01), long dashed line (a = 0.04), dashed-dotted line 
(a = 0.07), dashed line (a = 0.1), dotted line (a = 0.75) and continuous line (a = oo) (the 
exact solution of the Random Flight Model). 



Comment: Figure 8, First Author: Marcelo Marucho, JCP 



32 




FIG. 9. Characteristic function C (q) versus wave vector q for n = 20. 
Dashed-dotted-dotted line (a = 0) (the exact solution of rigid Model), 
dashed-dashed-dotted line (a = 0.01), long dashed line (a = 0.04), dashed-dotted line 
(a = 0.07), dashed line (a = 0.1), dotted line (a = 0.75) and continuous line (a = oo) (the 
exact solution of the Random Flight Model). 
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FIG. 10. Characteristic function C (q) versus wave vector q for a = 0.01. Dotted line 
(n — 6), dashed line (n = 10), dashed-dotted line (n = 20) and continuous line (n = 30). 



Comment: Figure 10, First Author: Marcelo Marucho, JCP 
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FIG. 11. Mean squared end-to-end distance (R 2 ) versus the parameter a for n = 8. The 
dotted line is our approximate solution and the dashed line is the exact solution of the KP 
model[? ]. 
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FIG. 12. Mean square end-to-end distance (R 2 ) (in logarithmic scale) versus the 
parameter a. Dashed-dotted line (n = 5), dashed line (n = 8) and dotted line (n 
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FIG. 13. ln{(R 2 )} versus the number of segments ln{n). Continuous line (a = 0.001) and 
dashed line (a = 0.75). 
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